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Abstract 


A set of relatively high resolution three-dimensional '3D) simulations were produced to 
investigate the triggering of moist convection by landscape generated mesoscale circu- 
lations. The local accumulated rainfall varied monotonically (linearly) with the size of 
mdn idual landscape patches, demonstrating the need to develop a trigger function that 
is sensitive to the size of individual patches. A new triggering function that includes the 
effect of landscape generated mesoscale circulations over patches of different sizes con- 
sists of a parcel s perturbation in vertical velocity (w 0 ). temperature (# 0 ). and moisture 
(q 0 ). Each variable in the triggering function was also sensitive to soil moisture gradi- 
ents. atmospheric initial conditions, and moist processes. The parcel s vertical velocity, 
temperature, and moisture perturbation were partitioned into mesoscale and turbulent 
components. Budget equations were derived for 6 0 and q 0 . Of manv terms in this set of 
budget equations, the turbulent, vertical flux of the mesoscale temperature and moisture 
contributed most to the triggering of moist convection through the impact, of these fluxes 
on the parcel s temperature and moisture profile. These fluxes needed to be parameter- 
ized to obtain 6 0 and q,.. The mesoscale vertical velocity also affected the profile of «\. 
We used similarity theory to parameterize these fluxes, as well as the parcel's mesoscale 


vertical velocity. 



1 Introduction 


Analysis of observational data and simulation results obtained over heterogeneous land 
surfaces can reveal both turbulent and mesoscale processes. For example, turbulent ed- 
dies can be superimposed upon mesoscale circulations that are generated by landscape 
discontinuities (Mahrt. 19S7. Balling. 19SS. Segal et at. 19S9: Smith et al.. 1992: Malirt 
et al., 1994: Chen and Avissar. 1994: Pielke et al.. 1997: Vidale et al.. 1997: Lynn et al.. 
199S). These mesoscale circulations generate sea-breeze-like fronts, which can be asso- 
ciated with relatively large vertical velocity, temperature, and moisture perturbations. 
S. udies with numerical models and some observational results have shown that these 
sea-breeze-like fronts produce both shallow and deep convection [e.g.. Chen and Avissar. 
1994: Cutrim et al. 1995: Avissar and Liu. 1996: Lynn et al.. 199$). 

In the past, parameterizations of atmospheric moist processes have been developed 
by dividing cumulus clouds into two groups, non-precipitating and precipitating (Frank. 
19S3). The former are shallow clouds that have a vertical depth of about 3 km. while 
precipitating clouds extend vertically into the middle- and upper-troposphere, and often 
form in very unstable atmospheres fZawadzki et a!.. 1981: Frank. 1983). Parameteriza- 
tions of shallow clouds usually assume that these clouds are triggered by Rayleigh-Benard 
convective instability (he., turbulent-scale processes (Wetzel and Boone. 1995) i. Param- 
eterizations of deep clouds, for the most part, relate the triggering of these clouds to 

large scale forcings, e.g.. the magnitude of the grid-scale moisture convergence (Rain and 
Fritsch, 1992). 

Rogers and Fritsch [1996) discuss a relatively new type of triggering function that 
includes the impact of landscape generated mesoscale circulations on the triggering vari- 
ables. They estimate the magnitude of the largest subgrid-scale vertical velocitv pertur- 
bation originating in each layer, and then calculate whether this perturbation is strong 
enough to overcome the total grid- resolvable negative inhibition between the source laver 
and the level of free convection. The work of Rogers and Fritsch (1996) is an important 
step in the process of developing a cumulus parameterization that is sensitive to the 
scale of the surface heterogeneity. However, they ignore subgrid-scale temperature and 
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moisture perturbations forced by tlie landscape heterogeneity 

Here, we use a high resolution cloud resolving model to simulate the atmospheric re- 
sponse to heating of landscape patches (described in Section 2 ). The purpose is multifold: 
i i To suggest a new approach to the triggering of landscape generated moist convection, 
and ii) to identify an appropriate set of triggering variables, (iii) Analyze the impact of 
landscape and atmospheric forcings on the triggering variables, and iv) show that the 
triggering variables relate well to the distribution of rainfall, vl Obtain mesoscale and 
turbulent components from the high resolution simulated data, and vi) calculate vari- 
ous flux terms in budget equations for the parcel's temperature and moisture. Then, to 
address the closure problem, an analysis of terms in budget equations shows which of 
the terms in these equations need to be parameterized. Finally, vii) Develop a parame- 
terization for these terms, as well as the mesoscale contribution to the parcel's vertical 

xelocity. These subjects are discussion in Section 3. A summary and conclusion is present 
in Section 4. 

2 Method 

2.1 Numerical Model 

A data set derived from observations that could be used for a study of the tvpe presented 
here does not exist. Such a study requires very high resolution data over a relatively 
large area (1000 km‘i. Fortunately, recent improvements in computer power have made 
possible a large number of high resolution model simulations. Some of these simulations 
have been applied to the study of convective processes over heterogeneous land surfaces 
(Nichols et ai. 1991: Chen and Avissar. 1994: Doran and Zhong. 1995: Fankauser it ai. 
1995: Avissar and Liu. 1996: Lynn et ai. 1998). Still, we note that observations like those 
with aircratt ' t.g.. houng. 1987. 19S8: Finkele et at.. 1995: Mahrt et ai. 1994: Seilers >t 
ai. 1995) ana Doppler radar (Houze. 1989) can provide the data required to verify model 
derived results and parameterizations. 

The moae: used in this stud}- is the Goddard Cumulus Ensemble Model 'GCE: Tao 
and Simpson. 1993). Briefly, the Rutledge and Hobbs (1984) scheme was used here to 
parameterize cloud microprocesses: the model uses a prognostic equation for the turbulent 
kinetic energy based on work by Deardorff (1975), Klemp and Wilhelmson , 1978). and 
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2.2 Experimental Domain 

The experimental domain had a 250x250 m* horizontal grid resolution with periodic 
lateral boundary conditions, and a stretched vertical coordinate (see Table 1 for more 
details). Each simulation was run for 12 hours with a time step of 5 seconds. The 
total domain size was 512 grid-elements in the west-to-east direction (or 12$ km) and 22 
grid-elements in the north-to-south direction (8 km). 

Two soundings on July 27 rA . 1991. one located on the east coast and one located on 
the west coast of Florida, were taken from CaPE (the Convection and Electrification Ex- 
periment I. They were averaged to obtain a mean sounding for an east-west cross section 
oxer the peninsula at 6 LST (Fig. 1). The sounding had a small initial con\-ective avail- 
able potential energy (CAPE) 2 of 740 k.J kg" 1 , but a relatively low lifting condensation 
lexe] pressure (an LCL of 1010 mb with a surface pressure of 101$ mb), loxv level of free 
convection (LFC: $39 mb), and high equilibrium level (EL: 190 mb). The vertical profile 
of the u-component of the wind is shown in Fig. lb (the v-component was set equal lo 
zero simply because we would expect it to haxe no important effect on ihe dexelopinu 
mesoscale circulations in the "truncated" north-to-south domain used here). The speed 
of the u wmd in the lower troposphere is small enough to allow for the development 
of landscape generated mesoscale circulations, but large enough to affect this develop- 
ment (see below). Upon moistening of the planetary boundary layer (PBL). the initial 
sounding was conducix e to the dex elopment of deep moist conx ection. 

We initialized the numerical model with a uniformly distributed homogeneous vec- 
etation consisting of a broad-leaf and coniferous forest. The stomatal resistance of the 
vegetation is a strong function of soil moisture in the root zone. For this reason, the 
soil moisture of sandy-clay-loam was chosen to "control" the surface distribution of land 
surface fluxes, leading to the generation of mesoscale circulations in our simulated do- 
mains. The process by which landscape circulations are generated by contrast imr patches 
of wet and cry ground has been described in detail by Pielke d a!. (1991 i. Axissar and 
•The CAPE was calculated using a mixed layer depth from the land surface to 500 m 
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riien (1993). and Lynn ft nl.. i 1995a I. We will not discuss this further and rite reader is 
reierred to these papers for further details. 

Seven land surface domains were used here to provide the surface boundarv conditions. 
In each, there were (different) alternating distributions of patches of dry and wet ground. 
In the dry ground, the volumetric soil moisture was chosen to be 10% of the fraction of 
the difference between field capacity and wilting point 3 . Here, this corresponds to a soil 
moisture. 8 si = 0.153. In the wet ground, this fraction was set at 90% of field capacity 
\0 S ; = 0.2S2). Superimposed on each grid element was a random value of volumetric soil 
wetness, r. which was —0.02 < r < 0.02. Land surface patches with 0 s! = 0.153 are 
referred to as "dry patches.' while land surface patches with 6 st — 0.2S2 are referred to 
as "wet patches." 


Four of the seven domains are shown in Fig. 2. and are labeled Domain 1 - Domain 4. 
In the main body of the text, we discuss in detail only results from Domain 1 and 
Domain 3. while in the Appendix we also refer to Domains 2 and 4. In Domain 1. there 
is a 64 km dry patch surrounded by two 32 km wet patches. Because the domain is 
periodic, there is actually a 64 km dry patch and 64 km wet patch in this domain. Note, 
the patches within the domains vary only in the west-to-east direction. In comparison, 
there were a number of differently sized dry patches in Domain 3. They are /left lo 
right) 7.5 km. 4 km. 16 km. 4 km. 8 km. 20 km. and 4 km in length. Domain 2 
contains patches of intermediate size, while Domain 4 contains very small patches. A 
summary of all domains used in this study is in Table 3. All domains are idealized 
representations of those that occur in situ. Still, the range of idealized domains chosen 
provides useful simulation results from which to draw conclusions about the potential 
impact of in situ landscape discontinuities on moist convection. In fact, the simulated 
impact of the landscape on the organized mesoscale flow was similar to that of Wane ft 
al. (1998). who used a statistical parameterization of the land surface heterogeneity. 

J V\ tiling point is the moisture content corresponding to matric potential of -15 bar. it is close to the 
volumetric water content at which plants can no-longer extract water from the soil. The field capacity is 
the volumetric water content corresponding to a balance that occurs in wet soil between diffusion (up. 
and gravitational drainage (down). 
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2.3 Numerical Experiments 

There were a number of experiments produced for this work. Seven simulations were 
pioduced using domains of various sizes with moist processes "turned on** i Exp. lw - 
Exp. < \\ ). These simulations were produced with the observed temperature and moisture 
sounding, and a background, westerly wind (constant with height) of 0.5 m s" 1 . Four 
more simulations were produced with the observed sounding, including the observed wind 
piofile shown in Fig. ib (Exp. Sw - Exp. llw). Two simulations were produced using 
Domains 1 and 3. but with moist processes “turned off and the constant background 
wind (Exp. 1. and Exp. 2). We produced two additional simulations with moist processes 
turned off and the observed wind sounding shown in Fig. lb (Exp. 3 and Exp. 4). By the 
phiase turned-on. we mean that clouds were allowed to form when model-grid elements 
reached supersaturation. By turned-off. we mean that no clouds were allowed to form in 
the simulated domain, and the atmosphere was allowed to super-saturate. Six additional 
simulations were produced with various initial boundary conditions. The reader can refer 
to Table 3 for more details. 

2.4 Budget Equations 

An\ model variable o can be separated into a resolved arid sub- resolvable component: 

o = o ~ o s ( 1 ) 

where the o represents the grid-scale, horizontal average at the large scale (i.t.. regional- 
or global-scale), and o s is the perturbation superimposed on it. We seek to develop 
a parameterization that can be used in regional- and global-scale atmospheric models. 
Thus, we assume that the flow component. o s ■ includes all motions on scales smaller than 
the synoptic flow. 

We find it useful to decompose o s into mesoscale. o', turbulent large eddv. o" . and 
turbulent small eddy. o'". components. Thus. 

o f = o' - o' 4- o'" , ■_> , 

Observations suggest that small, turbulent eddies have horizontal length scales of about 
oO m. while turbulent iarge eddies have horizontal length sale of about ] - 1.5 times 
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the boundary layer height \H: Hardy and Ottersten. 1969: Konrad. 1970: Kaimal ft a/.. 
1976: Caughev and Palmer. 1979: Taconet and Weill. 19S3i. The time-scale of these 
turbulent circulations is less than an hour. In comparison, mesoscale circulations have 
spatial scale of about 10 km - 200 km. and time-scale 50 min - 1 day [t.g.. Stull. 19SS: 
\ idale et a/.. 1997). Note, the decomposition of phi s is independent of the grid-scale of 
the cloud resolving model. 

We also assume that 

o = o -f o' + o" -f o'" 3 ) 

and 

(o) = (o + o' + o" + o'") = o (4) 

wheie the angle brackets (y indicate an average of the variable o l and its components i 
at the grid-scale of a regional- or global-scale atmospheric model. 

If the flow is one-dimensional, the variables, o. o'. o". and o'" can be formally defined 
from a discrete fourier transform [t.g.. Walker. 19SS) as: 
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H^re. A 0 <_i has a value that results in a filtering of both the large and small turbulent 
eddies from the data (i.e.. . the wavenumber Hf Ax - 1.5///AA' ('where A.r is the grid- 
spacing in the numerical model)). The field remaining after filtering is o + o'. from 
which o’ can be simply calculated. Likewise, the large eddies are found through the 
specification of an appropriate value for A>_, (he.. 50 m/A-ru In this case, the fields 
remaining would be o - o' 4- o" . Given both o and o', o" can also be calculated. 
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To produce budget equations for 6, and q s . we write a simplified equation for the 
perturbation o s . ignoring molecular diffusion (Stull ilfiSDii. 


do s - do s do do, 

-r <- >- fi/uT u 


dt 


dxj 


' d.r 


u j.s 


■ s dx, 


du^sOs 

dx, 


- •>* 


( 10 ; 


Ue then rewrite in flux form (we ignore horizontal advection of subgrid-scale quantities 
because we assume that circulations generated by patches are contained within each 
grid-element of the hosting model: fie., the regional- or global-scale model): 

do s dwOs , du' s o dw s o 5 


' dw\o. 


+ -Sa, 


( 11 ) 


dt d: d: ' dz \ dz , 

We then expand each term (but the first and last) by substituting for o s . We have 


dwo s die o' dil-o" du’o f ‘ 


d: dz ' d: dz 

du s o dw'o dw"o dw'"o 


dw s o s 

dz 


dz dz 1 dz dz 

dw'o’ du J o" dw'o'" dw"o' dw"o' 


dz 


dz 


dz 

dir'"o" dw'"o' dw'"o" dw"'o" 


dz 


dz 


( 12 , 


( 13 , 


(14l 


dz 


d: 


dz 


dz 


dw.o 


S^S 


d: 


dw'o' \ 

~ / 
dw"o " ' 
dz , 


dw'o" ] 


+ 


dz 

'dw"~'" 


dz 


/ dw'o'" \ 

\~dT / 

O'" \ / dw'"o ' ' 


dw"d 

~dT 


' dw'"o" \ / dw"' o’"' 


dz 


dz 


do) 


dz 


Finally, we can use the above equations to obtain (9, and q 0 . that is. a triggering 
parcels potential temperature and specific humidity. This is done by substituting 6 (and 

then q) for o. and then averaging the equation for 9, (and then q.) over the area of a 
parcel. Note. 


o s = o' + o" - o"" 


I lb j 


where the superscript^' indicates an average over the area of the parcel. 



3 Results 

3.1 A New Approach To the Triggering Problem 

In our opinion, two approaches have been suggested for triggering moist convection over 
heterogeneous land surface patches. The first approach assumes that the largest patch 
within a grid-scale domain produces the biggest sub-grid scale perturbations originating 
within each source la\er (Rogers and Fritsch. 1996). The second approach, as suggested 
by Lynn tt al. (19981. assumes that the size of the subgrid-scale perturbations should 
vary proportionally with the average size of the patches. In both, bigger patches should 
pioduce more rainfall than smaller patches because the perturbations over the former 
are larger than over the latter. 

Landscape dry patches heat the overlying atmosphere more than wet patches heat 
the overlying atmosphere. As a result, the atmospheric pressure drops over the drv patch 
arid mesoscale circulations can form on either side of the dry patch. These circulations 
mo\ed inward over the dry patch. In our simulations, mesoscale circulations occuring 
over individual dry patches produced rainfall over these patches. This rainfall occurred 
usually between 10 and 14 L.ST. 

We calculated the total accumulated rainfall from the simulations described in Table 3. 
We did so by integrating over the whole domain in time and space. A linear interpolation 
was then performed and the slope, y-intercept. and regression coefficient were obtained. 
The relationship between the domain averaged (accumulated ) rainfall and average patch 
size was not monotonic (Table 4). We obtained similar results even when we used the 
largest patch size within each domain, instead of the average patch size. Thus, neither the 

first or second approach would appear to well represent the rainfall that can be triggered 
by landscape patches. 

The reason that rainfall is not monotonic with average patch size is simple, and can 
be explained by two examples, shown in Fig. 3a (Exp. lwj and Fig. 3c (Exp. 4w f . After 
examining such figures, we can conclude that rainfall and its duration increase with 
increasing patch size 2 . However, isolated rain clusters can occur even over small patches. 

The forcing by the land surface depends upon the difference in heat flux between the wet and 
dry patches. Here, this difference at the time of rainfall formation was about 100 W m“' J , which can 
typically occur between different vegetation surfaces or in response t.o soil moisture differences if .a.. Sun 
and Maim. 1994) 
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ana ..remains (such as Domain 2> with small patches can have more patches iand hence 
more rainfall "clusters j than domains with large patches (such as Domain 1 i. Thus, the 
domain accumulated rainfall is a function of both patch size and patch number. 

h may be obvious to the reader that accumulated rainfall depends upon both patch 
size and patch number: but it might not be obvious how to proceed in the dev elopment of a 
triggering function. However, the results show a simple, linear relationship between patch 
size and accumulated rainfall over individual patches (Fig. 4j. In fact, the correlation 
coefficient was 0.99. 

Tae modeled relationship between rainfall and individual patch size can be explained 
b\ reterence to linear theory. In the absence of a background wind, linear theory has 
been used to simulate the coarse features of mesoscale circulations. Dalu ti a!. (1991) 
have mown that the intensity of the flow increases proportionally with increasing patch 
size, .or patches of size less than the local Rossby radius of deformation. 

Figures 3b and 3d show that the background wind profile can affect the development 
of rainfall over patches. A background wind increases turbulent dissipation, reducing, 
quite substantially. the generation of landscape generated rainfall over small patches. At 
the .'cme time, a synoptic wind from a warmer to a colder surface can strengthen the 
horizontal temperature gradient - provided, however, that the patch size is large enough 
and tr.e synoptic wind is less than a critical value - f.g.. 5 m s -1 (Pielke. Tims, 

tainfau oxer one side of the patch can increase in the presence of a background wind. 
a ~ shown in Fig. 3b. Thus, the distribution oxer individual patches, such as in Exp. Sw 
(Domain 1). can be quite asymmetric, suggesting the need to represent a multiple of cloud 
populations over such patches. Still, the correlation coefficient for the set of experiments 
(Exp. MV - Exp. 1 1 w ) was 0.95. 

Note, the model had cyclic boundary conditions. It seems clear that some of the 
ram occurred over the small dry patches relatively late in the simulation because of 
the boundary conditions. This was because some of the convective cells crossed the 
right boundary and reentered the domain on the left side. Had we used open boundarv 
condii.ons. the regression coefficient would likley have been even higher than obtained 
here. 

Tne nigh correlation coefficient in both cases suggests the validity of a third approach. 

A more realistic approach to the problem would be to develop a trigger function to be 
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applied to multiple, individual patches, whose associated variables would depend upon 
patch ?ize land atmospnenc background conditions). The domain averaged quantities 
could then be obtained from a cloud model as the integral over the domain of the regional- 
or global-scale atmospheric model of cloud related variables occurring over individual 
patches. 

3.2 Definition of Triggering Variables 

As noted above. Rogers and Fritsch (1996) recognized that sea-breeze-like fronts as- 
sociated with landscape generated mesoscale circulations can have parcels with strong 
(upward) vertical velocity. They empirically relate the size of the patches to the vertical 
velocity of the parcel, and use this velocity in their triggering function. 

To identify an appropriate set of triggering variables that could be used in a new 
triggering function, we examined modeled atmospheric fields at various times during the 
model simulations. These atmospheric fields were the modeled perturbation fields of hor- 
izontal wind (u,). vertical velocity (tr 5 ). potential temperature ($,). and specific humidity 
(?*)• 0ur analysis indicated, like that of Rogers and Fritsch (1996). that the most robust 
parcels occur in parcels along sea-breeze like fronts. However, it also showed that the 
temperature and moisture of these parcels can be different than the grid-scale. Thus, 
vertical velocity, temperature, and moisture affected the triggering of moist convection. 

For this reason, we propose a trigger function that includes the vertical profiles of 
velocity (uy). as well as potential temperature (i 9 0 ) and specific humidity (q 6 ). Here, the 
subscript zero means that the triggering variables refer to parcels moving along sea-breeze 
like fronts ( note, the parcel is an area average of the subgrid-scale perturbations). It seems 
evident that a trigger function that uses each of these variables has the advantage that it 
provides the information required to determine if triggering of moist convection should 
occur, without refering to an inhibition energy (Rogers and Fritsch. 1996) or empirical 
relationships between the grid-scale and subgrid-scale (Fritsch and Chappell. 1980) Still, 
out approach is similar to Fritsch and Chappell (1980) and Rogers and Fritsch (1996). 
in that we seek to relate subgrid-scale properties to the triggering of moist convection. 

This trigger function would be used as follows: i) The parcel s temperature and 
moisture would be used to calculate if saturation occurs on the subgrid-scale. In addition, 
the initial vertical velocity ol this parcel at the level of saturation would be used to more 
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realistically initialize a parcel s ascent, ii) Like Fritsch and Chappell (1980). The parcel 
uould ascend upward because it is buoyant: if it reaches its level of free convection, 
then convection would occur. Note, the triggering variables have a vertical profile that 
extends from the surface to the top of the PBL. Thus the parcel would grow into an 
environment that could be quite different (more moist and cool) than the srid-scale. 
The other schemes do not include a vertical dependence in the subgrid-scale variables. 

To obtain, tr-(-). &>(:). and qj : j for each dry patch, we first identified the location 
of the largest moisture perturbation along each front over each side of the patch. This 
became the center of our parcel for each frontal boundary. We then averaged u' s . ff s . and 
q s over surrounding grid-elements, to obtain parcels of square area (square areas were 
chosen here, rather than circular areas, out of convenience ). 

We examined the relationship between parcel size and the magnitude of the triggering 
variables. Figure o shows that there exists a strong sensitivity to the size of the parcel, 
but when the parcel size is either 10.6 or 1S.1 km 2 , its vertical profiles resemble those of 
the surrounding PBL (not shown). For our purposes, we chose the averaging size to be 
four grid elements, or a parcel of size 5.1 km 2 (9x9 grid-elements). This size is consistent 
v\ith the area of observed cloud bases, and those utilized in cumulus convection schemes. 

3.3 Impact of Initial Conditions on Triggering Variables 

We utilized a number of simulations with moist processes turned off. to study the depen- 
dence of the vertical profiles of the triggering variables on initial conditions. 

3.3.1 Sensitivity To Patch Size and Background Wind 

In Exps. lw and 4w. the rainfall occurred most intensely over bigger patches than small 
patches. To see why. we examined the triggering variables obtained with the same patch 
sizes, but with moist processes turned off (Exp. 1 and Exp. 2). When we averaged the 
triggering variables in Exp. 1 and Exp. 2 over this time period, we noted that, for the 
most part, the triggering variables over patches of different sizes had similar magnitude 
(Fig. 6). 5et. differences were obtained in rainfall over the different patches. Clearlv. the 
vertical depth of the triggering variables had an impact on rainfall. This vertical depth 
reflects the relative robustness of the mesoscale circulations over each patch. 
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In Exp. Sw and 11 u . the background wina i the observed sounding i had strong 
enect on rainfall. To see why. we examined results from Exp. 3 and 4. which have the 
'-sine patch sizes as Exp. Pw and llw. but moist processes were turned off. Figure 7 
shows the vertical profiles of the triggering variables along the downwind side of each dry 
parch (also between 10 - 14 LST). We show these profiles, rather than the profiles on 
the upwind side of the patches because rainfall was most intense along this side of the 
various patches. In this case, both the magnitude and vertical depth of the triggering 
variables depended upon patch size. 

A background wind increases the turbulent dissipation of the atmospheric gradient in 
temperature between the dry and wet patches. Over small patches, this dissipation was 
sufficient to reduce the magnitude of the triggering variables, even over the downwind 
side of each patch. Thus, in Exp. llw. the peak values of rainfall over the smaller patches 
were much less than in Fig. 4w. 

However, as noted above, a synoptic wind from a warmer to a colder surface can. if 
the dry patch is large enough, strengthen the horizontal temperature gradient. This can 
enhance the frontal forcing at the downwind patch boundaries. This affect is present in 
Exp. 3. and was reflected in the rainfall distribution of Exp. $w. 

To emphasize this point. Figure S provides a comparison of triggering variable results 
obtained on the upwind and downwind side of the dry patch in Exp. 3. It shows that 
the background had a detrimental effect on the triggering variables on the upwind side 
of the patch (the westward (left) side). This is because a synoptic wind direction from a 
colder to a warmer surface weakens the horizontal temperature gradient. Thus, the set 
of triggering variables obtained over a individual patch can depend upon the direction 
and magnitude of the background wind. We will show that it is useful to include such 
information in a parameterization of moist convection. 

3.3.2 Impact of Surface Forcing, Atmospheric Stability. Specific Humidity, 
and Latitude 

In Exp. 5. we examined the impact of the soil moisture gradients upon the triggering 
variables (Fig. 9i. The changes in soil moisture of the dry and wet patches affected 
the horizontal contrast in surface sensible heat fluxes, and subsequent vertical, turbulent 
transfer of heat and moisture (not shown). This led to relatively less robust mesoscaie 
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circulations than in Exp. 1 (the control experiment u and a decrease in the magnitude 
of die triggering variables. These variables decreased roughly proport ionallv with a 
dec lease in the standard deviation of the surface sensible heat fluxes (not shown i. Lynn 

( (1995b) also found a similar relationship between the mesoscale heat fluxes and the 
standard deviation of the surface sensible heat flux. 

In Exp. 6. we increased the atmospheric stability of the sounding. The new sounding 
was typical of a summertime day with strong subsidence occurring within a strong high 
pressure system. This impacted both the surface fluxes and atmospheric turbulent fluxes, 
reducing the depth of the vertical mixing when compared to Exp. 1. The change in 
stability strong]}' effected the vertical depth of the triggering variables. 

Moreover, the magnitude of w 0 obtained in Exp. 6 was much less than obtained in 
Exp. 1. The variable w 0 is high!}' sensitive to the turbulent convection, which depends 
strongly on the stability. However, the vertical convergence of the heat fluxes determines 
the vertical structure of B, and q Q , Here, these fluxes converged over a more shallow laver 
than in Exp. 1. Thus. 0 O and q 0 were affected less by the change in stability than ay. 

The percentage change in the initial profile of the domain average specific humidity. 
9* uas simply proportional to the change in q 0 (Exp. < j. Moreover, the change in q led 
to a decrease in the sensible heat fluxes from both dry and wet patches, which had a 
negative impact on the development of mesoscale circulations. Thus, the change in q 
aEo affected the magnitude of the vertical profiles of uy and 6L; each was smaller than 
in Exp. 1. 

The change in latitude had relatively little effect on the triggering variables (Exps. 8 
and 9). The patches were too small to allow for a significant impact of the Coriolis force 
on the triggering variables. 

3.4 Relationship between Triggering Variables and Moist Pro- 
cesses 

Above, we showed that the vertical depth and magnitude of the triggering variables 
corresponded to the intensity of rainfall over patches of different sizes. We did so l.v 
comparing profiles of the triggering variables from simulations without moist processes, 
and surmised that the results obtained in those simulations would apply to simulations 
with moist processes. Wt. it is useful to look in more detail at the relationship between 
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raintall and the triggering variables, using results from simulations with moist processes. 
We can provide further support for the previous conclusions, as well as emphasize addi- 
tional particulars of this relationship. To limit the length of the discussion, we show oniv 
vertical profiles of q 0 . 

The growth and decay in the vertical profiles of q., corresponded quite well with 
the growth and dissipation of rainfall in each experiment (Fig. 10). Similar results were 
obtained ior w 0 and 6 0 (not shown ). In addition, the time-scale of the triggering variables 
obtained over Patch #2 and #5 in Exp. 4w were much less than Exp. 2. Similarly, rainfall 
grew and decayed more quickly over these patches than over larger patches. 

In addition, the time scales over Patch #2 and ~r-5 were short enough to lead to a 
growth of q 0 . its dissipation, and then regrowth of q 0 (Fig. 10c. d t over these patches. In 
contrast, the growth and dissipation of the triggering variables over the larger patches 
took longer than over the smaller patches. Similarly, the rainfall occurred over these 
patches for a longer time period than over the smaller patches. 

Figure 11 shows that the triggering variable q c grew in size more quickly than its 
counterpart on the upwind side of the dry patch (Exp. Sw). The differences in q , obtained 
on the upwind and dowjn wind side of the dry patch corresponded with differences in 
rainfall shown in Fig. 3b on each side of this patch. 

Prior to the formation of moist convection, the trigger function variables are associ- 
ated with mesoscale circulations that are forced solely by the land surface patches. After 
the formation of moist convection, however, these circulations can be extensive! v modified 
by moist convection, t.g.. by cloud shading of the surface, and evaporative cooling along 
the sea-breeze-fronts (Lynn et al. 1998). A comparison of Figures 10 and 11 with Figs 6 
and i show the impact of these processes on the triggering variables. Moreover, turbulent 
processes occurring within clouds can transfer heat and moisture upwards into middle 
troposphere. These proceses also affect the vertical profiles of the triggering variables. 

Cloud shading also weakens the forcing on the developing mesoscale circulations. 
However, evaporative cooling counteracts, in part, the impact of cloud shading on the 
developing circulations (Lynn it al.. 1998). It also accelerates the movement of the 
fronts inward towards the center of the dry patch. Thus, mesoscale circulations affected 
by moist processes have a shorter time-scale than those unaffected by moist processes. 

A more thorough discussion of this issues is beyond the scope of this paper. 
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3.5 Budget Equations 

We analyzed the budget equation? using the data set described in the Appendix. The 
Appendix describes how the modeled data were separated into mesoscale \ o'} and large- 
eddy (o") perturbations. We did not obtain o'" from that modeled data because the 
resolution of this data was much bigger than the horizontal length scale of these eddies. 
Note, the model's parameterization for turbulence represents small eddy processes, and. 
of course, a proportion of the unresolved exchange due to larsre eddies. Converselv. lanje 
eddy fluxes calculated with the filtered fields might include contributions from small 
eddies. 

Because we are unable to obtain the small eddy perturbations, we have a simplified 
budget equation for o 0 . 

dwo s dwo' dwo" 

~ t — ! x — ! 1 7 S 

dz dz dz 

dw s o du-'o du-"o 

~dT == ~dT + ~dT tI -' 

dlL' s O s duz'o' du-'o" dw"o' du-"o" dw'"o'" 

T — ^ ~ 5 ! o 1 3 i ^ 1 1 9 ) 

dz dz az dz dz dz 

du^oA 

dz ] 

where the overline indicates a Reynold's average over the time- and -spatial scale of 
the small eddies. Note, the Reynolds assumption implies that the fluxes wo”’. w"'o. 
do'". w"'o'. and w''o'" = 0. A test of these equalities is not formally presented here, 
and is left for interested parties. However, we note that the parameterized turbulent 
fluxes were large only in the lowest part of the PBL. This suggests that the cross terms 
described above are indeed small above the lower PBL. and that the calculated large-eddy 
perturbations do indeed represent the large-eddy variables that might occur m stiu. 

Fig. 12 shows the vertical profiles of various terms in the budget equation, obtained 
in Exp. 3. We chose Exp. 3 because it used the observed sounding (with the u wind), 
which is. perhaps, more typical than a light-wind case. In the budget equations, there 
are a number of fluxes for which parameterizations do not yet exist. Vet. there is a need 
to develop a parameterization for only u-"i9 ,P and w' 'q' y . These terms are the vertical 
transport of the mesoscale temperature and moisture by the turbulent wind. New pa- 
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rameterizations for other terms nor previously described in the literature were found not 
to be required, t.q.. \c'9' f . because these terms were relatively small. 

These results can be explained through an analysis of Figure 13. which shows the 
vertical profiles of the mesoscale and large-eddy perturbations obtained in Exp. 3. It 
is for a "triggering' parcel moving upwind (east-to-west ). with an upwind moving >ea- 
breeze-like front. The data were averaged over the time period 10 - Id LST. Quite simplv. 
w" >> u J . while 9" « 9'. and q" << q' . 

Figure. 14 shows the vertical profiles of the fluxes in the budget equations, but for 
those obtained in Exp. Sw. A careful examination of this figure suggests the followins: 
i l Both ir"9' H and ic"q' p again contribute most significantly to the triggering of moist 
convection, ii) The vertical structure of these curves remained relatively unchanged, 
thus, we can develop a single parameterization for both w"9 and w"q'' that should 
work prior to and during the formation of clouds, iii) Both u 77 ¥ p and t TUf*' force 
the vertical transport of heat and moisture within the cloud, iv) These terms, which 
represent the vertical transport of the t urbulent potential temperature and moisture, can 
be parameterized using conventional cumulus cloud parameterizations. Table 5 shows a 
listing of terms that occur in the budget equations. It also shows which of these terms 
are required to obtain the vertical profiles of 9 0 and q,,. 

The results obtained above were found to apply to other experiments produced as 
part of this work. Therefore, we believe the conclusions reached to be valid for initial 
conditions other than those examined here. Note, we do not show in any of our figures 
the fluxes ( w'"9 { w"'q"') . or xc l "9'" P . w"'q"' H but these contributed significantly to the 
vertical profiles of 9 0 and q 0 within the surface layer and lowest 100 or so meters of the 
PBL. Nor. do we show the vertical profiles of Se or S q . Each should be obtained as a 
residual of the cloud parameterization. 

3.6 Parameterization 

We made a thorough examination of model results obtained in the previous experiments. 
W'e concluded that the method suggested by Lynn it at. i 1995b) should be used to obtain 
parameterizations for the needed terms. The reason being that the parameterized fluxes 
iand u J i have a time-scale different than the turbulent fluxes. Thus, to realistically 
describe the evolution of the fluxes requires a prognostic equation for them, rather than 
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a diagnostic relationship between them and the surlace land atmospheric) forcing. We 
applied Buckingham Pi theory to obtain the dimensionless numbers for the parameteri- 
zation (Stull. 19SS). 

A Chebyshev polynomial has been used to describe the vertical structure of the fol- 
lowing dimensionless numbers: 
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where a k . and 3 k are empirical parameters deduced from the experimental (numerical) 
data. set. and r is linearly related to the time of day /. 

The coefficients are determined by solving numerically Eq. 24 in two steps, using a 
time-splitting technique (e.c.. Pielke. 1954). In the first step, intermediate coefficients 
c" are calculated using an explicit scheme. 
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and in the second step, the c k are calculated using an analytical solution: 

c;-i " - Ar) = c k (r + 


(26) 


To represent the residual impact of the heat fluxes on the PBL. D x as an e-folding time 
of three hours. Note, we have a set of three c k . cjf . cf e ‘. and cf q '. as well as a set of 
three a k . which are of. a k " e ' . 


and q{. 5 . We refer to onlv one 3 k . where 
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The dimensionless number ZT consists of: 
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while D- is: 


D : = f(R 0 . L) . 2'* . 

The Appendix describes the form of the parameterization in more detail. Table 6 has a 
listing and description of all variables used in the parameterization, and Table 7 has a 
listing of the empirical coefficients a K . 

We used model output to obtain the complete set of dimensionless numbers. The 
sensitivity of the dimensionless numbers, including the fluxes and mesoscale vertical 
velocity, to changes in surface and atmospheric variables was similar to that obtained for 
the mesoscale heat fluxes by Lynn et al.. 1995b. The reader is referred to their work for 
more details. Instead, we only present Fig. 17 and IS. which shows the vertical profiles 
of D^jrpp. D^yrpp. and D—p from the model output and the parameterization for each. 

These figures demonstrate the ability of the parameterization to represent the sen- 
sitivity of the model results. However, the magnitude of the fluxes obtained were less 
than simulated by the model. This is because of the non-linearity of the fluxes ( <r" varies 
chaotically in time). Still, the parameterization captured the most important features of 
the observed data set: i) the dependence upon patch size, ii) background wind, iii) sta- 
bility. and iv) specific humidity. Other cases from the experimental data are not shown, 
but were independent!}' tested to insure that the parameterization for each variable is 
robust. 

4 Summary and Conclusion 

A set of relatively high resolution three-dimensional (3D) simulations were produced to 
investigate the triggering of moist convection over heterogeneous land surface domains. 
This moist convection was triggered by mesoscale circulations generated by the landscape 
heterogeneity. We found that a monotonic (linear) relationship exists between the local 
accumulated rainfall over individual patches and the size of these patches, but not the 
domain accumulated rainfall and domain averaged patch size. Thus, we suggest that cu- 
mulus parameterizations and their trigger functions for heterogeneous landscapes should 
be applied over multiple, individual patches within the domain, rather than to a single, 
patch of average size. 

An appropriate set of triggering variables for this triggering function consists of the 
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vertical profiles of vertical velocity, temperature. and moisture, for parcels moving a lorn* 
sea-breeze like fronts associated wit!: landscape generated mesoscale circulations. We 
suggest a simple trigger function using these variables, which uses the principles described 
by Fritsch and Chappell (19S0) and Rogers and Fritsch (1996). 

The triggering function variables were found to be sensitive to the surface contrasts 
in sensible heat flux, atmospheric stability, initial profile of atmospheric moisture, and 
moist (cloud) processes. They were relatively insensitive to latitude, bur showed that 
background wind could impact the magnitude of the triggering variables. a.$ well as 
create important differences between triggering variables over individual patches. We 
noted that the evolution of the triggering variables prior to and during the development 
of moist convection described quite well the local distribution of rainfall over patches. 

We derived budget equations for 9 and q 0 . These equations contain contributions 
from mesoscale and turbulent fluxes, as well terms such as w"9' y and ir"q ,p \ which repre- 
sent the turbulent vertical flux of the mesoscale temperature and moisture perturbation, 
respectively). We used a Fourier transform to filter the data, and obtain a distribution 
of mesoscale and turbulent perturbations. Using this new data set. we then calculated 
the various terms in the budget equations, and determined that u J, 6 ,v and ic"q ,? need 
to be parameterized to close each respective equation. Quite interesting] v. parameteri- 
zations developed for cases without ciouds should also apply to cases with clouds, since 
the vertical structure of these terms remained relatively unchanged by moist processes. 

It might be interesting to the reader to note that the grid-scale mesoscale fluxes are 
also insensitive to phase transitions. Therefore, the parameterizations developed bv Lvnn 
et al. (1995b) and Zeng and Pielke *1995) should also apply as given for both drv and 
moist regimes. 

As noted, a linear relationship was obtained between rainfall and patch size. Perhaps, 
if moist convection had occurred later in the simulation, the relationship found would 
have a more exponential shape (increasing non- linearly upward with increasing patch 
size). This is because non-linear advective affects become more important relative to 
linear processes as mesoscale circulations develop and move inward over the domain. 
Moreover, we would expect that the maximum rainfall would occur for patch sizes equal 
to the local radius of deformation i(.c.. Lynn ti al.. 199$). The sizes of patches simulated 
here were all less than the local radius. 
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We used similarity theory to parameterize w"0' p and w"q' p and tr'. The important 
variables that were incorporated into a parameterization are the patch size (and local 
radius ot deformation), gradient in the surface sensible heat fluxes, the background wind, 
planetary boundary layer height, and planetary boundary layer atmospheric potential 
temperature and specific humidity. These variables were used to develop empirical rela- 
tionships between them, the fluxes, and the mesoscale vertical velocity. 

Finally, the development and dispensation of data sets such as FIFE and LBA can 
provide additional evaluation and refinement of the parameterizations suggested in this 
work. 


5 Appendix 

5.1 Filtering the Modeled Data 

To analyze the budget equations, we needed to obtain from the modeled data the 
mesoscale and turbulent fields. Observations suggest that turbulent eddies hat e horizon- 
tal length scales of about 1-1.5 times the boundary layer height (Hardy and Ottersten. 
1969: Konrad. 1970: Kaimal ft ai. 1976: Caughey and Palmer. 1979: Taconet and Weill. 
19S3I. It is a simple task to calculate the boundary layer height, and then obtain the 
appropriate wavenumber for filtering the data. Note. Table S contains a summarv of 
definitions used in this section and elsewhere in this paper. 

We calculated the boundary layer height across the domain and used the maximum 
boundary layer height to obtain a filtering wavenumber for the mesoscale fields. We 
then examined the vertical profiles of the triggering variables obtained after assuming 
that turbulent eddies have wavelength from 1.5 to 3.5 (in increments of 0.5) times the 
boundary layer height i H ). The distribution of the atmospheric fields were most sensitiv e 
to filtering in the range of values 1.5 H - 2.5 H. However, the distribution of the filtered 
fields did not change very much when we used a wavelength of 3.0 H or 3.5 H . This implies 
that when we used, for example, a wavelength of 3.0 times the boundary layer height, 
that this value enabled a filtering of the large-eddies from the modeled data, leaving the 
mesoscale fields. 

To support this supposition, we did a spectral analysis ; Fig. 15). using the domains 
shown in Fig. 2. This analysis revealed that the dominant spectral energy (in the u- 
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wind field ■ corresponds to the wavenumber of the surface forcing. For example, peaks 
in the power spectrum occurred, for example, for wavenumber 4 (f.g.. Patch = 3: a 
•’>2 km dry-wet patch couplet) and wavenumber 7 (Patch =5: a IS km couplet). Most 
importantly, a simulation of very small patches (=s 0.2-5 m*: Domain 4i. produced an 
energy spectrum with wavelengths smaller than about 10 km. The energy associated 
with wave:engths greater than 10 km. was mesoscale-kinetic energy ( Avissar and Chen. 
109-3). It corresponds to the set of mesoscale perturbations derived from the initial data 
set. In contrast, turbulent kinetic energy occurred in association with wavelengths less 
than 10 km. This energy corresponds with the large eddy perturbations. Note, the 
scale-separation between the two regimes depends upon patch size. 

Based upon the prior analysis, we adopt a wavelength of 3.0 H to filter the data. The 
value of 3.0// obtained from our simulations differs from the range \H - 1.5 H suggested 
by observations. This :s because the model simulations used a relatively coarse resolution 
to resolve the turbulent fields. As noted, the smallest eddies that the model could resolve 
with a minimum of two grid-elements corresponded to a wavelength of 500 m. This is 
much larger than the typically observed wavelength for the small eddies. 

Figure 16 shows an example of filtered atmospheric fields obtained in Exp. 3. The 
mesoscale Horizontal wind (u'j. potential temperature (O'), and specific humidity fields </ 
corresponded very well with observations of sea-breezes and sea-breeze-like circulation- 
produced with mesoscale models ' f.g.. Finkele ti al.. 1994: Lynn ft al.. 1995a). However, 
the model simulated vertical wind field was "contaminated" by the buoyancy forcing as- 
sociated with turbulent motions. After filtering in space (and time) the modeled vertical 
velocity field, we could not obtain a meaningful mesoscale vertical velocity field. For this 
reason, we do not show this field in Fia. 16. 

Insteac. we note that the hydrostatic assumption applies when the ratio of the vertical 
to the horizontal length scale of the circulation is equal to or less than one i. Pieike. 
19S4). Thus, we believe it justified to assume that the derived mesoscale fields were 
hydrostatic. Moreover, the mesoscale vertical velocity field obtained from the mesoscale 
horizontal wind (using the continuity equation) corresponded well with results obtained 
from mesoscale models and obtained indirectly through the use of observations. We >how 
this derived field in Fig. 16b. and obtain the vertical velocity in al! experiments from the 
continuity equation. 
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As described previously. Fig. 13 shows vertical profiles of the mesoscale and turbu- 
lent perturbations obtained for the triggering parcels in Exp. 3. The relatively robust 
mesoscale perturbations occurred in this experiment because of the contrast in turbu- 
lent transport between the dry and wet patches, and horizontal and vertical advert ion 
associated with the mesoscale circulations. In contrast, turbulent, large eddies formed 
because of buoyancy forcing. Near the ground, the difference between these eddies and 
their surroundings was relatively small. Hence, the magnitude of their perturbations near 
the ground was relatively small compared to the mesoscale perturbations. In the upper 
PBL large eddy thermals are relatively cool and moist when compared to entraining air 
from rhe capping inversion from above the mixed layer. This entrainment generated more 
significant turbulent eddy perturbations in the upper PBL than occurred in the lower 
PBL. 

For the purposes of the analysis, we also obtained the turbulent, large eddy fields. 
For brevity, these are not shown. However. o n = o s — o' . 


5.2 Equations for Parameterization 


The dimensionless number. D\ is denned as follows: 
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The dimensionless variable A xr"'6"' is used to represent the relative forcing bv the land 
surface. The bigger it is. the stronger this forcing. It increases with the growth in the 
difference in sensible heat flux between patches, but decreases when U > L' 0 . 

The function is defined for two types of parcels: those moving with the 

prevailing wind (/((■’.£'“)) and those moving against the prevailing wind {f\V.L'f > . To 
obtain f(l .[ “). we first define the ratio of {’ to where V is the background wind 
and L" u is a constant = 1.0 m ? -i . 

Then, for parcels moving with the prevailing wind: 
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To obtain f(l\ L f ). we also define the ratio of (.' to where U* varies as a function 
of parch length. To obtain this functional relationship, we made additional simulation*; 
with Domain 1 (Table 3). In one case, a sea- breeze-like front was nearly stationary over 
the downwind side of the dry patch. Here, the background wind (in the PBL) was about 
4 m s -1 . Most interestingly, the vertical velocity of the vertical wind perturbation at 
this stationary front (prior to the convergence of the fronts' was most robust of any of 
the simulations. Bechtold tt al. (1991) discuss the reasons why stationary sea breezes 
produce the most robust vertical velocities. To simplify the application of their results 
to the parameterization, we simply assume that U 0 — 6.2-5E-5 s -1 L (= 4 m s _1 for 
Domain 1). The dependence of l 7 on patch size allows for the impact of turbulent 
mixing on the horizontal gradient in temperature. Thus, we have: 
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For parcels moving against the prevailing wind: 
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For D' 2 - we define the variable .r as: 
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and 0 < x 0 < 1 Note, the variable R 0 is a local radius of deformation (see Lynn d 
al.. 1998). The frictional dissipation variable in the equation for R 0 was determined to 
be three-hours from an examination of the mesoscale kinetic energy obtained with the 
cumulus ensemble model. 

For all cj 
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d A = do 

dn = d 2 

(37) 
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Figure 1: Atmospheric sounding obtained from the Convective and Precipitation Elec- 
trification Experiment. (CaPE). (a) Skew-T diagram of temperature and dewpoint: 
(b) profile of the initial, u component of the background wind. 

Figure 2: Distribution of land surface domains in some of the numerical experiments. 

Figure 3: Cross-sectional plots of rainfall versus time obtained in (a) Exp. Iw. ( b i 
Exp. 8w, (c). Exp. 4w, (d) and Exp. llw. 

Figure 4: Total accumulated rainfall versus patch size for the dry patches listed in Table 4. 
Exps. Iw - 7w. 

Figure 5: The vertical profiles of the triggering variables obtained after averaging over 
square areas, centered at the location of the maximum moisture perturbation on the 
downwind side of the dr}' patch in Exp. 3. The label of each curve is the size of the 
square parcels with units of km 2 . The parcels represented were located along the 
sea-breeze like fronts over the dry patch. 

Figure 6: Vertical profile of a\ . 6 ;V . and q 0 obtained in Exp. I and Exp. 2 (for labeled 
patches in Fig. 2. Domain 3). The data presented represent the average of profiles 
obtained from 10 - 14 LST. 

Figure 7: Same as Fig. 6. but for Exp. 3 and Exp. 4 (for labeled patches), obtained on 
the downwind side of each dry patch. 

Figure S: Same as Fig. 6. but for the triggering variables obtained on the upwind side 
and downwind of the dry patch in Exp. 3. 

Figure 9: Same as Fig. 6. but for Exp. 1 and the sensitivity experiments described in 
Table 3 \ the number on each curve corresponds to the number of the experiment in 
Table 3 i. 
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Figure 10: \ertical profile of q G obtained in Exp. lvv and 4w averaged over four hourly 
periods. These hourly periods are: fa.) 10 - 11 LST. and ibi 11 12 LST. (ri 12 - 

13 LST. and (di 13 - 14 LST. The arrow indirat.es the approximate hieght of the 
cloud base. 

Figure 11: Same as Fig. 12. but for the triggering variables on the upwind side and 
downwind of the dry patch in Exp. 8w. 

Figure 12: Vertical profile of (a) heat and (b) moisture fluxes obtained in Exp. 3. The 
vertical convergence of these fluxes (not shown) contributes to the time rate of change 
of and q Q . The labels on the profiles in (a) correspond as follows: (1) u J 9 ,v . (2) 
ii J, 6' fP . (3) u J 9” r . (4 i w”9 fP . (5) {w’9')* and (6) {ic n 9"). The labels on the profiles 
in (b) are as follows: ( 1 ) w l q ,p . (2) ir"q" p . (3) (4) u Jf q ,t . 15) n ,f q f ;. and lb) 

\w n q n ). The o v indicates an average over the area of the parcel. Moist processes were 
turned off in this experiment. 

Figure 13: Mesoscale (meso) and turbulent (turb) perturbations obtained in Exp. 3 for 
the triggering parcel on the downwind side of the dry patch. The data shown were 
i he average of these perturbations obtained between 10 and 14 LST. 

Figure 14: Same as Fig. 12. but for Exp. Sw. Moist processes were turned on in this 
experiment. 

Figure 15: Power spectrum analysis of the u component obtained in experiments using 
the domains indicated (see Fig. 2). Simulations were produced using the observed 
sounding, with a light background wind (0.5 m s^ 1 ). The data were averaged over 
12 hours of simulation. 

Figure 16: Vertical west-to-east cross section of mesoscale perturbation fields obtained 
at 13:00 LST in Exp. 3. Note, the derived mesoscale perturbations were independent 
of y . or the north-to-south direction. The dry patch was located from 32 to 06 km. 
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Fig. 17: Vertical profiles of the triggering variables obtained from model simuia; ion result^ 
between 12 14 LST * left side) and 16 - 18 LST i right side) from experiments 

indicated by number on curve. 

Fig. 18: Vertical profiles of the triggering variables obtained from the parameterization 
between 12 - 14 LST (left side) and 16 - IS LST (right side) from experiments 
indicated bv number on curve. 


31 



List of Tables 


Table 1: Model input parameters used for numerical simulations (Numbers in parent 
refer to simulations with moist processes "turned on.” 

Table 2: Land characteristics used for 1 lie numerical simulations. 

Table 3a: Description of model experiments with moist processes turned on. The r 
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Table 3b: Description of model experiments with moist processes turned off. The r 
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. the correlation coefficient between the average patch size and rainfall was 0.77. 
while the correlation coefficient between the biggest patch size and rainfall was 0.70 
\ in these two sets, the range of the rainfall was relatively small). 

Table 5: Budget equations terms. An A" in the appropriate column indicates that the 
corresonding term needs to be included in its particular budget equation. A yj in- 
dicates that an appropriate parameterization needs to be developed for that term; 
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an average over the area of the triggering parcel. o Ci — o fV -f ~o^ P . The variable o' 
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Table 1: Model input parameters used for numerical simulations (Numbers in paren- 
theses refer to simulations with moist processes "turned on.*’ 


Condition 

Yclut 

Day of the Year 

•July 27 

Latitude 

•j v 

Initialization time 

6 a. m. 

Integration time step 

o s 

Simulation length 

12 hr 

Height of the atmosphere 

10 km 20 km) 

Number of vertical grid elements 

30 50) 

Vertical Grid Resolution (Stretched) 

20 - 500 m 

Lateral Boundary Conditions 

Periodic 

Horizontal Grid Resolution (Fixed) 

250> 250 m 2 
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Table 2: Land characteristics used for the numerical simulations. 


Land Ch a ra ctn -is / i c 

\ alvc 

Surface Roughness (m) 

i.u 

Surface Albedo 

0.20 

Surface Emissivity 

1 .00 

Soil Texture 

sandy-clay-loam 

Soil Depth (m) 

10 

Root Zone Depth (m l 

1.60 

Porosity (min 3 mm” 1 ) 

0.404 

Field Capacity (min' 3 nun" 3 ) 

0.29S 

Wilting Point (mm 3 mm" 3 ) 

0.137 

Slope of Retention Curve 

6.77 

Saturated Matric Potential (cm) 

-13.49 

Saturated Hydraulic Conductivity (cm -1 ) 

4.45 xlO -4 

Vegetation 

Tall. Broadleaf 
and Needleleaf Trees 

Leaf Area Index 

2.46 

Minimum Stomata.] Resistance (s/m) 

1.00 

Maximum Stomata! Resistance (s/m) 

10.0 


Plant Critical Water Potential (m) 


-200 




Table 3a: Description of model experiments will) moist processes turned on. I he r 


component of the wind was set equal to 0. 


Xann Patch Si~cs (L. hnj Special (' audition 


Exp. lw 04 (Domain 1 ) 


Exp. 2w 

40. 32 

Exp. 3w 

16. 24. 24.5 (Domain 2) 

Exp. 4w 

7.5. 4. 16. 4. 8. 20. 4 (Domain 3) 

Exp. 5vv 

3 < L < 8 

Exp. 6w 

2 < L < 0 

Exp. 7w 

0.25 < L < 1 (Domain 4) 

Exp. 8w 

64 

Exp. 9w 

16. 24. 24.5 

Exp. 10w 

40. 32 

Exp. llw 

7.5. 4. 16. 4. 8. 20. 4 


Background u wind a constant. 0.5 m s 
Background u wind a constant. 0.5 m s 
Background u wind a constant. 0.5 in s 
Background u wind a constant. 0.5 m s' 
Background u wind a constant. 0.5 m s' 
Background u wind a constant. 0.5 m s 
Background u wind a constant. 0.5 in s' 
Observed v background wind 
Observed u background wind 
01)served u background wind 
Observed ?/ background wind 



Table 3b: Description of model experiments with moist processes turned off. 1 lit' r 


component of the wind was set equal to 0. 


A dine 

Patch Sizes (L . km > 

Spt rial ( 'and Hum 

Exp. 1 

64 (Domain 1 ) 

Background u wind 0.5 m s -1 

Exp. 2 

7.5.4. 16. 4. 8. 20 (Domain 3) 

Background a wind 0.5 m s" 1 

Exp. 3 

64 

Observed i u ) background wind 

Exp. 4 

7.5.4. 16. 4. 8. 20 

Observed i u) background wind 

Exp. o 

64 

6 S { — 0.177 (dry patches). 0 s l = 0.257 (wet patches! 
Background v wind 0.5 m 

Exp. 6 

64 

Increased stability by 3.5 1\ 1000 m” 1 
Background v wind 0.5 m s~ ] 

Exp. 7 

64 

Decreased q to 50% of observed 
Background a wind 0.5 m s _1 

Exp. 8 

64 

10 u latitude 

Background v wind 0.5 m s -1 

Exp. 9 

64 

50° latitude 

Background v wind 0.5 m s~ ] 

Exp. 10 

64 

Observed ( u ) equals one-half background wind 

Exp. 1 1 

64 

Obsen*ed (w) equals twice background wind 

Exp. 12 

16, 24. 24.5 (Domain 2) 

Background u wind a constant. 0.5 m s~ ] 

Exp. 13 

0.25 < L < 1 (Domain 4) 

Background a wind a constant. 0.5 m s -1 
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Table 4: Rainfall versus average patch size and largest patch size. Exps. lvv - Exp. 7\v 
were each produced using a light background wind. Exps. Sw - lhv were produced using 
the observed (west-to-east) background wind. The correlation coefficient for average 
patch size versus rainfall for Exps. lw - Exp. 7w w'as 0.32, while the correlation coefficient 
between the biggest patch size and rainfall was 0.45. For Exps. Sw - llw, . the correlation 
coefficient between the average patch size and rainfall was 0.77, while the correlation 
coefficient between the biggest patch size and rainfall was 0.79 (in these two sets, the 
range of the rainfall was relatively small). 


Name 

Average Patch Size (km) 

Biggest Patch Size 

Accumulated Rainfall (mm) 

Exp. lw 

64 

64 

1.35 

Exp. 2w 

36 

40 

1.2S 

Exp. 3w 

21.5 

24.5 

1.31 

Exp. 4w 

9 

20 

1.48 

Exp. 5w 

^ 4.5 

si 10 

1.33 

Exp. 6w 

ss 3.0 

si 6 

1.23 

Exp. 7w 

Si 0.5 

Si 1.0 

0.98 

Exp. 8w 

64 

64 

1.18 

Exp. 9w 

36 

40 

1.13 

Exp. lOw 

21.5 

24.5 

1.14 

Exp. llw 

9 

20 

1.14 
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Table 5: Budget equations terms. An X in the appropriate column indicates that the 
corresonding term needs to be included in its particular budget equation. A \J indicates 
that an appropriate parameterization needs to be developed for that term; otherwise, 
parameterizations exist for that term in the literature. The ^ indicates an average over 
the area of the triggering parcel. d> 0 = 4? P + . The variable 4>' represents the mesoscale 

field obtained after a filtering of the high resolution data. The variable 4>" represents the 
large eddy turbulent field obtained as the residual of the filtered and original fields. 
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Term 


Include Develop Parameterization 


< w'O ') 


(w"8") 

A' 

(w'"0"') 

X 

w'8' P 



X 

w ff 6 fP 

X 

wW* 


w"'6" ,P 

X 

W) 


W) 

A 

W) 

X 

w f q tV 



X 

w"q' V 

X 

w l q i,P 


u?"V" p 

X 
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Table 6: Variables required for Parameterization. 


Variable 

Description 

Units 

L 

Length Scale of Patch 

m 

Ro 

Local Radius of Deformation 

km 

H 

Planetary Boundary Layer Height (for each dry patch) 

m 

'3 

<1 

Gradient of Surface Sensible Heat Fluxes 

I< m s -1 

e 

Grid-Scale (Mean) Planetary Boundary Layer Temperature 

K 

Q 

Grid-Scale Specific Humidity 


u 

Background Wind 

m s' 1 

u u 

A function of patch size for parcels on upwind 



side of dry patch 

m s" 1 

u d 

0 

For parcels on upwind side of dry patch 

111 s- 1 

tau 

time 

s 
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Table 7a: Constants required for the vP 9 . 


Variable 

Value 

Gi 

SO. 3 

q 2 

-76.6 

alphas 

-73.0 

a 4 

75.8 

Os 

20.2 


Table 7b: Constants required for the w /f 9 /P . 

Variable Value Units 



0.55 

a 2 

-0.41 

alphas 

0.20 

a 4 

0.53 

<*5 

-0.42 


Table 7c: Constants required for the w"q fP . 

Variable Value Units 


Oi 

32.1 

<*2 

20. S 

alphas 

-31.1 

a 4 

-26.5 

O 5 

21.9 
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Table 8: Miscellaneous Definitions of Variables ( H is the planetary boundary layer 
height). 

Variable Definition 

0 Grid-scale average over domain of regional model 

<j>* Mesoscale perturbation (horizontal-scale > than 3 H) 

0" Large eddy perturbation (50 m < horizontal-scale < 3 H) 

<f) fn Small eddy perturbation (horizontal-scale < 50 m) 

4> s Subgrid-scale perturbations obtained with GCE 

0o Parcel value of 0; 0 S averaged over square area of parcel 

0^ Average of 0 over square area of parcel 
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U) U) 






Domain 1 



32 64 96 128 

Distance (km) 


Domain 3 


2 3 5 6 



32 64 96 128 
Distance (km) 


Domain 2 



Distance (km) 


Domain 4 



Distance (km) 



0.16 0.18 0 20 0.22 0.24 0.26 0 28 0.30 

Soil Moisture (v/v) 


0.10 0.12 0.14 





















































- 1.0 - 0.5 0.0 0 . 

Potential Tempe: 





o 





sqrt(Power) (m/s) 





HEIGHT (km) HEIGHT (km) 


I k 




-6 -4 -2 -1 - 0.5 - 0.2 0.2 0.5 1 2 

SCALE 



















- 0.25 0.00 0.25 0 . 5 < 

w' p (m s -1 ) 



- 1.0 - 0.5 0.0 0.5 

WW P (k m s- 1 ) 



0.5 










